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We prove variationally that at weak coupling in one, two, 
and three dimensions there exist correlated electron-phonon 
states below the approximate ground states characteristi- 
cally found by adiabatic polaron theory. Besides differing 
non-trivially in quantitative aspects such as the value of the 
ground state energy, these improved ground states are found 
to differ significantly in qualitative aspects such as correla- 
tion structure and scaling behavior. These differences are 
sufficiently severe as to require a reevalutation of the phys- 
ical meaning attached to such widely used terms as "large 
polaron" and "self-trapping transition". 

PACS numbers: 71.38.+i, 71.15.-m, 71.35.Aa, 72.90.+y 

The theory of polarons in condensed matter is a ven- 
erable subject that has drawn the attention of many cre- 
ative minds applying an impressive array of theoretical 
techniques for more than half of this century. In more re- 
cent days, there has been a surge of activity in this area 
driven in part by a revitalized interest in the potential 
role of polarons in superconductivity and by some pivotal 
advances in computational techniques. 

Although there is much now that can be done to reveal 
the secrets polarons have held close for many years, many 
of the concepts shaping the perspective from which these 
are viewed are decades old. It should not be surprising 
that some updating might be required. 

One such concept is the dichotomy of "free" vs. "self- 
trapped" states and the intimately related issues of the 
dependence of polaron structure on lattice dimensional- 
ity. This concept arises from adiabatic theory [|l] pd| , 
which we use rather broadly to include a number of 
semi-classical approaches and uncontrolled approxima- 
tions that may not explicitly invoke the adiabatic ap- 
proximation in an obvious way, but which are essentially 
indistinguishable in their result. According to adiabatic 
theory, the minimum energy state of an electron-phonon 
system of any dimensionality at sufficiently strong cou- 
pling is a compact, spatially-localized state. With de- 
creasing electron-phonon coupling strength, however, the 
energy of this branch of localized solutions is found to rise 
until, for D > 2, it penetrates the free-electron energy 
band at some finite coupling strength characteristic of the 



particular lattice structure and type of electron-phonon 
interaction. On the weak-coupling side of this crossover, 
the minimum energy states are found to be the free elec- 
tron states; hence the notion of a self-trapping transi- 
tion arises as a crossover in the global ground state from 
"free" to "self-trapped" character. In 1-D, the minimum- 
energy solution at any value of electron-phonon coupling 
is found to be a localized state; hence, the notion that all 
1-D electron-phonon states are " self-trapped" . 

Although appealing in many ways, very widely ac- 
cepted, and correct in certain respects, almost every "in- 
teresting" element of this adiabatic characterization of 
polaron structure and self-trapping is here proven to be 
break down in the weak-coupling regime, even in 1-D. 
This breakdown of adiabatic theory can be demonstrated 
with a simple and now decades old variational method at- 
tributable to Merrifield @^3). Although the Merrifield 
method is not the method of choice for wide-ranging, 
high precision calculation, its use for our present purpose 
lends some transparency to the demonstration since it is 
quite closely related to some very typical adiabatic ap- 
proaches while simultaneously comparing favorably with 
more accurate methods Q . 

We use the Holstein Hamiltonian [ p^|Jl6| 

H = - J 4i°* + nw b rfi" 

< m , n > n 

- 9^ a n a n( b n + b rt) , (1) 
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in which at creates a single electron in the rigid-lattice 
Wannier state at site ft, and &t creates a quantum of 
vibrational energy Hlo in the Einstein oscillator at site 
ft. All sums run over the entire D-dimensional lattice of 
edge-length N; the restriction < m,ft > limits rfi and ft 
to be nearest neighbors along the each crystal axis, and 
J is the electron transfer integral between such neighbor- 
ing sites. We note that throughout this paper, "coupling 
strength" refers to the dimensionless local coupling pa- 
rameter <?, and not the relative measure A = g 2 fiuj/2J 
common in adiabatic approaches. 

In one very characteristic approach, the ground state of 
the electron-phonon system is approximated by a Pekar 
form 



1 



m = \a) 



(2) 



in which |a) is a one-electron state and |/3) represents 
the (multi-phonon) state of the lattice. The lattice state 
1/3) is sometimes subsumed into a classical treatment of 
the lattice; however, since the lowest-energy lattice states 
consistent with arbitrary classical configurations are co- 
herent states, a very typical choice for the components of 
the Pekar form are 



exp 
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Applying such a state to the model Hamiltonian and then 
varying the parameters {a,/?} to minimize the energy 
yields the adiabatic character outlined above. That is, in 
all of 1-D and at sufficiently strong coupling in 2-D and 
3-D, a balance can be struck between the spread of the 
electronic component (a) and the focussing tendencies of 
the lattice component to yield non-trivially localized 
ground states, while at weak coupling in 2-D and 3-D, 
the free electron ground state is energetically favored. 

Rather than account for electronic spreading within 
the Pekar form, the Merrifield method uses a delocalized 
state that strictly satisfies the Bloch condition regardless 
of coupling regime 
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This state can be obtained from (Q) by first localizing the 
electronic component (setting a n = S n o) and then form- 
ing the appropriate superposition of displaced replicas 
(k — for the global ground state) . 

As with the Pekar form, the Merrifield state is ap- 
plied to the system Hamiltonian and the parameters 
{/3} are varied to minimize the energy. The resulting 
self-consistency equations are then solved numerically. 
Our calculations were performed in nearly-real time on 
a single-processor Sun Microsystems Sparc 10 desktop 
workstation. Although for the parameter values consid- 
ered here convergence to "bulk" values is accomplished 
on modest-sized lattices, the data presented in this pa- 
per were all obtained on isotropic cubic lattices of edge 
length N = 64 in order to obviate any potential concerns 
relating to boundary effects; thus, in 1-D, 2-D, and 3-D, 
N D = 64, 4096, and 262, 144 respectively. 

Although the Merrifield method is not accurate at 
intermediate and strong coupling (its distorts the self- 
trapping transition and converges toward the first order 
of strong-coupling perturbation theory, missing the im- 
portant second-order correction p7|), it has the virtue 
of providing a straightforward variational confirmation 



of weak-coupling perturbation theory. This is signifi- 
cant because the nonlinear nature of the stability argu- 
ments supporting the adiabatic characterization of po- 
laron structure naturally raise some concerns over the 
reliability of weak-coupling perturbation theory that are 
hereby put to rest. 

As the comparisons in Figure [l] show, the Merrifield 
method provides variational proof that in the weak- 
coupling regime: i) the free electron states are not the 
lowest energy states in any dimension, 2) the "large po- 
laron" states of 1-D adiabatic theory are not the low- 
est energy states in 1-D, and 3) the correlated electron- 
phonon states that do minimize the energy are more 
consistent with weak-coupling perturbation theory than 
with adiabatic theory. This consistency between weak- 
coupling perturbation theory and the Merrifield method 
is not limited to the ground state energy, but extends to 
the polaron effective mass (LSI, kinetic energy |l9| , and 
the spatial structure of the electron-phonon correlations 
as well J2(l] ; examples of the latter are shown in Figure |^. 
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FIG. 1. Ground state energy in units of Two in 1-D, 2-D, and 
3-D for J = hu). Chain-dotted line: Bottom of the free elec- 
tron band at —2D J. Dotted line: Strong-coupling asymptote, 
E(0) ~ — g 2 . Dashed line: Weak-coupling perturbation the- 
ory. Bold solid line: Merrifield method. Faint solid line: Adi- 
abatic theory on a discrete 1-D lattice; data kindly provided 
by G. Kalosakas [ilp^ . Points: Quantum Monte Carlo, data 
kindly provided by P. E. Kornilovitch |23|J 
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We are thus forced to consider the crucial question of 
just how to characterize the correlated electron-phonon 
ground states of the weak-coupling regime. For histor- 
ical reasons at the heart of this paper, there has been 
considerable reluctance to characterize these states as 
"large polarons", though in every dimension these are 
the weak-coupling complement to the strong-coupling 
states that are unambiguously and universally charac- 
terized as " small polarons" . The principal reason for 
this reluctance appears to be that "large polaron" has 
come to be identified with the archetypical finding of adi- 
abatic theory in 1-D that the minimum-energy states at 
weak-coupling are broad, pulse-shaped structures char- 
acterized by a locking relation between the electron and 
phonon amplitudes (e.g., (3 n = g\a n \ 2 when using (|l|) - 
(^)) and by certain scaling properties in the continuum 
limit. The Pekar form (H)-(|) m particular applied to (|l|) 
at sufficiently weak coupling (permitting the continuum 
limit) yields 



a{x) 
(3(x) 



1/2 



sech(Ax//) , 



sech 2 (A;c/Z) , 



(G) 
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in which A = g 2 Ticu/2J and I is the lattice constant. 

Yet, the results displayed in Figure [l] prove that "large 
polaron" ground states as here described do not exist 
in 1-D for the same reasons that "free electron" ground 
states do not exist in 2-D and 3-D; viz., the 1-D "large po- 
larons" of adiabatic theory lie at higher energy than the 
correlated electron-phonon states indicated with great 
mutual consistency by weak-coupling perturbation the- 
ory, band-theoretic variational methods, quantum Monte 
Carlo, and others Q| . 

This disagreement is not a casual matter of complex 
correlations being described with fewer or greater pa- 
rameters, or exponential vs. gaussian tails, or differing 
numbers of retained orders that ultimately are of little 
consequence when handled with appropriate perspective 
and care. Rather, the adiabatic notion of the "large po- 
laron" in 1-D differs from actual weak-coupling polaron 
structure in the most essential of polaron properties, the 
polaron size and in the scaling of this size with system 
parameters. For example, consider the correlation func- 
tion 
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that can be viewed as measuring the shape of the polaron 
lattice distortion around the instantaneous position of 
the electron, or, essentially equivalently in view of the 
strongly local character of the electron-phonon coupling, 
as an image of the electron density one associates with a 
localized polaron. 

In 1-D, using the collected results following from the 
Pekar form, the spatial variance of this correlation func- 
tion is straightforwardly found to be 
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which we note diverges strongly at weak coupling. On 
the other hand, the result that follows analytically from 
weak coupling perturbation theory in any dimension pc| ] 
is 
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which is independent of g and scales differently with 
J/hu than does the adiabatic result. This result is con- 
firmed at weak-coupling by the Merrifield data in 1-D, 
2-D, and 3-D (see Figure g), and by data from the more 
accurate Global-Local method in 1-D (higher D pending) . 
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FIG. 2. Correlation decay along the x axis in 1-D, 2-D, and 
3-D, as measured by Cl 23 ' /Cq D \ Symbols indicate the results 
of weak-coupling perturbation theory. Segmented curves in- 
dicate the results of the Merrifield method for g — 1/2. Solid 
lines: 3-D. Chain-dotted lines 2-D. Dashed lines: 1-D. Bold 
lines, circles: J = 1. Faint lines, diamonds: J = 8. 

For the cases illustrated in Figure ^, the "large po- 
larons" of adiabatic theory would be considerably wider; 
based on the square root of the ratio of the variances, 
about 5 times wider in the J = 1 case and 14 times 
wider in the J = 8 case. 

These results reflect a weak-coupling polaron structure 
that is more rigid than is expected by adiabatic theory 
in 1-D, but more correlated with phonons than is ex- 
pected by adiabatic theory in 2-D and 3-D; spatially, 
weak-coupling polarons are much more compact than is 
expected by adiabatic theory in 2-D and 3-D, but may 
be larger or smaller than the "large polarons" expected 
by adiabatic theory in 1-D. 

Rather than stumble over jargon associated with dis- 
tinctions now evaporated, it would appear most sensible 
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to regard the qualitatively-similar correlated electron- 
phonon ground states of the weak-coupling regime as 
large polarons regardless of lattice dimensionality. More- 
over, since we have demonstrated that "free" states do 
not exist as the minimum-energy states at weak cou- 
pling in any dimension, the notion of self-trapping as a 
change in the ground state from "free" to "self-trapped" 
character does not express the physical content of the 
phenomenon. The self-trapping transition nonetheless 
exists as a clear feature of the polaron landscape. In 
more operational and model-independent terms it is the 
more-or-less rapid transition from states characteristic of 
the weak-coupling regime to states characteristic of the 
strong-coupling regime as reflected in observable polaron 
properties; i.e., a transition from large polaron structure 
to small polaron structure. Again contrary to adiabatic 
theory, which asserts the absence of a self-trapping tran- 
sition 1-D because of the absence of free-electron ground 
states, this notion of a self-trapping transition does exist 
in 1-D as well as in 2-D and 3-D. 

It is not the case, of course, that adiabatic theory is 
everywhere invalid. It appears from multiple considera- 
tions, however, that part of what occurs in the course of 
the self-trapping transition is a dematerialzation of the 
locking relation between electron and phonon coordinates 
that is essentially universal in adiabatic approaches [BSl . 
On the strong-coupling side of the transition, this locking 
is endemic and semi-classical treatments appear to en- 
joy respectable agreement with fully-quantum mechan- 
ical strong-coupling perturbation theory; on the weak 
coupling side, however, while still ultimately local in 
character as reflected in Figure |[ electron-phonon corre- 
lations are increasingly sensitive to the discrete quantiza- 
tion of phonon energies and momentum-space structure. 

These finer details are not well-captured quantitatively 
by the Merrifield method, but a very accurate and de- 
tailed description is available in the Global-Local varia- 
tional method that follows from direct, sequential gen- 
eralization of the Merrifield method. We have used the 
Global-Local method to analyze 1-D polaron structure in 
detail throughout the polaron parameter space E3 , 56[ , 
and higher-dimensional calculations by this method are 
in progress ]2^] . 

ACKNOWLEDGEMENT 

The authors gratefully acknowledge P. E. Kornilovitch 
and G. Kalosakas for providing numerical data used in 
Figure |l|. This work was supported in part by the 
U.S. Department of Energy under Grant No. DE-FG03- 
86ER13606. 



[1] Y. Toyozawa, Prog. Theor. Phys. 26, 29 (1961). 

[2] Y. Toyozawa, in Polarons and Excitons, edited by C. G. 

Kuper and G. D. Whitfield (Plenum, New York, 1963), 

p. 211. 

[3] D. Emin, Adv. Phys. 22, 57 (1973). 

[4] A. Sumi and Y. Toyozawa, J. Phys. Soc. Japan 35, 137 
(1973). 

[5] D. Emin and T. Holstein, Phys. Rev. Lett. 36, 323 
(1976). 

[6] M. Ueta et al, Excitonic Processes in Solids (Springer- 

Verlag, Berlin, 1986). 
[7] V. V. Kabanov and O. Y. Mashtakov, Phys. Rev. B 47, 

6060 (1993). 

[8] Y. Toyozawa and Y. Shinozuka, J. Phys. Soc. Jap. 48, 
472 (1980). 

[9] E. A. Silinsh and V. Capek, Organic Molecular Crys- 
tals: Interaction, Localization and Transport Phenomena 
(AlP-Press, New York, 1994). 

[10] K. S. Song and R. T. Williams, Self-Trapped Excitons 
(Springer Verlag, Berlin, 1996). 

[11] Y. Toyozawa, Pure & Appl. Chem. 69, 1171 (1997). 

[12] R. E. Merrifield, J. Chem. Phys. 40, 445 (1964). 

[13] Y. Zhao, D. W. Brown, and K. Lindenberg, J. Chem. 
Phys. 106, 5622 (1997). 

[14] A. H. Romero, D. W. Brown, and K. Lindenberg, to ap- 
pear in J. Chem. Phys. (1998). 

[15] T. Holstein, Ann. Phys. (N.Y.) 8, 325 (1959). 

[16] T. Holstein, Ann. Phys. (N.Y.) 8, 343 (1959). 

[17] A. S. Alexandrov and S. N. Mott, Polarons & Bipolarons 
(World Scientific, London, 1995). 

[18] A. H. Romero, D. W. Brown, and K. Lindenberg, cond 
mat, 9809025 (1998). 

[19] A. H. Romero, D. W. Brown, and K. Lindenberg, cond 
mat 9806031 (1998). 

[20] A. H. Romero, D. W. Brown, and K. Lindenberg, cond 
mat 9808348 (1998). 

[21] G. Kalosakas, S. Aubry, and G. P. Tsironis, Phys. Rev. 
B 58, 3094 (1998). 

[22] G. Kalosakas, private communication (1997). 

[23] P. E. Kornilovitch, cond-mat 9808155 (1998). 

[24] P. E. Kornilovitch, private communication (1998). 

[25] D. W. Brown and K. Lindenberg, Physica D 113, 267 
(1998). 

[26] D. W. Brown, K. Lindenberg, and Y. Zhao, J. Chem. 

Phys. 107, 3179 (1997). 
[27] A. H. Romero, D. W. Brown, and K. Lindenberg, in 

preparation (1998). 



4 



